Spin swap gate in the presence of qubit inhomogeneity in a double quantum dot 
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£Sj ' We study theoretically the effects of qubit inhomogeneity on the quantum logic gate of qubit swap, 

which is an integral part of the operations of a quantum computer. Our focus here is to construct a 
robust pulse sequence for swap operation in the simultaneous presence of Zeeman inhomogeneity for 
quantum dot trapped electron spins and the finite-time ramp-up of exchange coupling in a double 
^— i , dot. We first present a geometric explanation of spin swap operation, mapping the two-qubit 

\^ • operation onto a single-qubit rotation. We then show that in this geometric picture a square- 

' pulse-sequence can be easily designed to perform swap in the presence of Zeeman inhomogeneity. 

Finally, we investigate how finite ramp-up times for the exchange coupling J negatively affect the 
performance of the swap gate sequence, and show how to correct the problems numerically. 

d ' PACS numbers: 03.67.Lx 

The study of quantum information processing has attracted a great deal of attention in recent years because of the 
potential advantages provided by quantum mechanical principles such as superposition and entanglement . Among 
the many proposed quantum computer (QC) architectures, solid-state schemes such as those based on electron [2|, y] 
and nuclear spins Q and those based on superconducting circuits @ are widely regarded as the best candidates in 
providing scalable systems. However, solid-state architectures also have their own shortcomings. For example, in 
supeconducting-circuit-based QC schemes, offset charge noise due to jumps of trapped charges has been an important 
limiting factor [Hill; In semiconductors such as silicon, complexities in their band structures can potentially lead 
to significant difficulties in qubit manipulations 0. The often small energy scales in solids lead to slower coherent 
operations and harder initializations. Furthermore, in many solid-state QC schemes it is almost impossible to have 
completely identical qubits. As potential and problems coexist in proposed solid-state QC architectures, theoretical 
explorations are sorely needed to provide effective forewarnings and guidelines. In particular, it is important to analyze 
CO . theoretically various possible sources of errors (arising, for example, from imperfection, inhomogeneity, decoherence, 
nonadiabaticity, and in general from deviations from ideal architectures assumed in the QC proposals which must 
invariably be present in real solid-state systems) in these solid-state QC schemes. For example, In the spin-based 
quantum dot quantum computer (QDQC) architecture pi IToL ITlL Il2| . where trapped electron spins are the quantum 
bits (qubits), quantum dots provide the tags and the environment for the individual qubits. Each quantum dot is 
generally slightly different in size, geometry, confinement potential depth, g factors, etc. Some of these differences can 
CO ' be accounted for straightforwardly by system calibration, while others, such as inhomogeneity in the electron spin 
, Zeeman splitting, have to be treated more carefully. 

In this paper we study how to overcome the problems arising from qubit inhomogeneity and gate imperfections 
in the spin-based QDQC, particularly focusing on the swap operation. It is important to emphasize here that in 
realistic QC architectures, two-qubit operations include not only the entangling operations such as controlled-NOT, 
1 , \ but also auxilliary operations such as swap, which are crucial components for the effective manipulation of a QC 
(the function of swap in QDQC is to move non-neighboring qubits together and apart through a quantum dot array 
during entangling operations). Furthermore, controlled-NOT operation in QDQC is built upon the square-root-of- 
swap operations. Therefore, imperfections that affect swap will also affect controlled-NOT in general. For instance, 
previously we have shown that swap operation jp cannot be done precisely in one step if there is inhomogeneity in 
the electron Zeeman splitting [l^, because the inhomogeneity breaks the symmetry of the two-spin system. Here 
' we first present a geometric interpretation of two-spin swap operation in terms of single spin rotation, and design a 
square pulse sequence to perform swap in the presence of Zeeman inhomogeneity. We then discuss the effects of qubit 
- - ■ inhomogeneity and imperfect control of the exchange coupling on the swap gate sequence, and present its numerically 
corrected version. 

Under the condition that the double quantum dot low-energy dynamics can be described by Heisenberg exchange 
Hamiltonian [lTl[T^| . the two-electron spin Hamiltonian can be written as a sum of the exchange term and the Zeeman 
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splittings: 

H = .m 1 -s 2 + 11 s lz + l2 s 2z , (i) 

where J gives the strength of exchange coupling and is a function of quantum dot size, geometry, and confinement. 71 
and 72 are the Zeeman splittings of the two spins and are functions of local g factors and magnetic fields. If we express 
this Hamiltonian on the two-spin basis | TT), | Tl), | IT), and | J.J.), and let S = 71 — 72 be the Zeeman inhomogeneity 
and 7 = (71 + 72)/2 the average Zeeman coupling, we obtain 
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H= n 7 * n > (2) 



Notice that within this Hamiltonian the two unpolarized states | Tl) and | J,|) are decoupled from the two polarized 
states I TT) and | |j), thus their dynamics can be described separately. 

To understand how swap works, we first write down the spin states before and after a swap: 

a i "\ o, ( 

0*2 



£J = ai/3i|tt) + «2/3 2 |U) + ai^|U) + a 2 /3iUT) 

a rr , a rr , a 1 /3 2 + a 2 f3 1 , ai/3 2 - "2^1 
= aipiTf + a 2 p 2 li H = J 



V2 

ft ) ^ ( «2 ) = aiAI TT) + a2/?2 ' U) + a2AI U) + ai/32 ' U> 

= ai/^iJ T + a 2 P2Ji H -y| Jo o, 

where T T = | TT), T; = | ||), and T = (| Tl) + I it))/V2 are the triplet states, while S = (| U) - I iT))/v^ is the 
two-spin singlet state. From the above expressions, swap is achieved by switching the coefficients of the unpolarized 
I Tl) and I |T) states, or equivalently, changing the coefficient of the singlet component by a ir phase shift relative to 
the triplet states. The phase shift can be easily obtained by a Heisenberg exchange Hamiltonian H = JSi ■ S 2 , since 
singlet and triplet states are the eigenstates of the exchange Hamiltonian (split by 4J) so that its only effect is to 
introduce dynamical phases to each basis state. 

If a uniform magnetic field is present (71 =72 =7), the triplet states are split by the Zeeman coupling 7. When 
the exchange Hamiltonian is applied to the two-spin system together with the uniform magnetic field, swap can still 
be performed, with an additional phase: instead of an exact swap, now the final states take on the form 

a 2 ) ® U J - { ) ^ [ a 2 e~^ ) ' (3) 

In other words, different spin states acquire different phases depending on their Zeeman energies |13|. Thus the 
additional phases here come only from the polarized states. 

When there is Zeeman inhomogeneity between the two quantum dots, the spin dynamics is more complicated 
since the singlet and unpolarized triplet states are coupled by the inhomogeneous Zeeman splitting. However, as we 
mentioned before, within Hamiltonian Q the dynamics of the polarized and unpolarized two-spin states are separated. 
With the polarized states still eigenstates of Hamiltonian (J), we need to focus on only the unpolarized states | Tl) 
and I IT)- Since these states are not coupled to the polarized states, the two-state subspace they span can be treated 
as an effective spin- 1/2 system with an effective Hamiltonian [l5| 

He = ( j J gj= J(J x + = a ■ ( 4 ) 



where J n = %/ J 2 + 8 2 and n — (J, 0, 5)/J n . Within the picture of this effective spin-1/2 system, there is a simple 
geometric explanation to the two-spin swap operation. Hamiltonian Q is a rotation on the Bloch sphere of the 
effective spin-1/2 system around the axis given by h, with the rotation angle being determined by the duration and 
strength of this Hamiltonian. In the absence of inhomogeneity (5 = 0), the rotation is around the x axis. Thus a 7r 
rotation for this effective spin-1/2 system would correspond preceisely to a "swap" in the original two-qubit system: 
I Tl) - I IT) and I |T) - I Tl)- 
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In the presence of Zeeman inhomogeneity, the rotational axis h deviates away from the x axis. Now starting from 
the north pole of the Bloch sphere, the state will not be able to reach the south pole by one rotation around a 
fixed axis (corresponding to a square pulse of exchange coupling in the presence of inhomogeneity). Thus, Zeeman 
inhomogeneity makes the exact swap impossible by a single pulse of J (with a fixed sign) |13|. However, swap can still 
be performed if the Zeeman inhomogeneity is known (for example, if it is due to engineered g-factor). Essentially, we 
can adjust the rotational axis by changing the magnitude of the exchange coupling J. For example, we can switch on 
a square pulse of exchange with magnitude J for a tt rotation, then turn off the exchange and let the system undergo 
a tt rotation around the z axis (driven by the inhomogeneous Zeeman splitting) , then switch on a second square pulse 
of exchange with magnitude J' = <5 2 / J for a tt rotation. If we define R(h, 6) to be a ^-rotation around the n direction, 
the pulse sequence for swap would be: 

[/swap = R(n, 7r) R(z, 7t) R(n ', 7r) , (5) 

where the rotational axis are n = (J,0,S)/VJ 2 + S 2 , z = (0,0,1), and n' = (J', 0, S)/V J' 2 + 5 2 . The end result of 
this pulse sequence would be a swap for the two spins: | jj,) «-» | J.|). Figure n snows a schematic comparison of 
the square pulse sequence in the absence and presence of Zeeman inhomogeneity. If in the above pulse sequence the 
exchange couplings J or J' is unphysically large, we can always use smaller exchange couplings instead, but with 
a longer pulse sequence. For example, if J is the maximal exchange while J' is even bigger (which implies a large 
Zeeman inhomogeneity 5), we can keep performing R{n, 7r) R(z, tt) until at last a R(h' , tt) with J' < J can be used to 
complete the sequence (J' is determined uniquely by J and 5): 

U swap = R(h, tt) R(z, tt) R(h, tt) R(z , tt) • ■ ■ R(n', tt) . (6) 

Take a simple example of J/5 = tan(7r/8) = v2 — 1 with J being the largest possible physical value, the pulse sequence 
for swap is 

C^swap = R{h 1 tt) R(z, tt) R(n, tt) R(z, tt) i?(n, tt) R(z, tt) i?(n, tt) , (7) 

which invokes the exchange coupling 4 times. This simply shows that exchange becomes less efficient in performing 
swap operation if Zeeman inhomogeneity is large. An alternative interaction might have to be used instead (such as 
optically assisted spin flip). 

Using the geometric picture for swap, it is apparent that the presence of inhomogeneous Zeeman splitting introduces 
additional complexities to swap operation. For instance, in the absence of Zeeman inhomogeneity, it does not matter 
whether the exchange coupling J is switched on suddenly or gradually as the rotational axis is always fixed along 
x; while in the presence of Zeeman inhomogeneity, swap operation cannot be done for two spins (or, flip cannot 
be achieved for the effective spin-1/2) in one shot. The introduction of a finite pulse rise/fall time aggravates this 
problem as the orientation of the rotational axis of the effective spin-1/2 system becomes time-dependent (a gradually 
switched-on exchange coupling means that the rotational axis is time-dependent for the switch-on period, which 
immediately leads to incomplete and/or imprecise rotations). To quantify this difference we perform a numerical 
calculation and show that the shift in rotational axis reduces the efficiency of exchange coupling in performing swap 
operations. We then numerically search for the appropriate exchange coupling to perform swap operation under these 
non-ideal conditions. 

In the following we use the square pulse case (sharp rise and fall of the exchange coupling) as a benchmark to 
measure the degradation of the rotation by the exchange pulses with finite rise/fall times and introduce an effective 
infidelity Aer z = cr 2 | opt — o~q, where <7o is the smallest o~ z reached by a square exchange pulse after a tt rotation starting 
from the north pole of the Bloch sphere, while cr 2 | op t is the smallest er z reached by an exchange pulse with finite r r 
and optimized pulse duration (not necessarily a tt pulse, which is ill-defined for a time-dependent axis anyway). Fig. [2 
gives several trajectories for optimized rotation on the Bloch sphere projected onto the xy plane. Each trajectory 
represents a case with a chosen pulse rise/fall time and an optimized pulse duration so that the final destination point 
has the smallest (cr z ) component (closest to the (cr z ) — —1 south pole point. Recall that a rotation from the north 
pole to the south pole and vice versa corresponds to an exact swap: | fj.) «-> | IT))- All the curves share a common 
exchange coupling J — 0.2meV and an inhomogeneity of 5 — O.lmeV. One interesting feature here is that when the 
pulse rise/fall time is sufficiently long (5ps in this case), the system has to undergo approximately one and a half full 
rotation (instead of one half full rotation, or the 7r-pulse) in order to reach the smallest (o~ z ). In addition, notice that 
these optimal rotations are generally not 7r-rotations anymore. 

In Figures [3J and [SI we show how the efficiency of an exchange pulse in performing rotations varies with the 
exchange coupling J, Zeeman inhomogeneity 6, and the pulse rise/fall time r r . In these figures, the filled and unfilled 
symbols correspond to two different pulse shapes: sinusoidal and linear rise/fall, respectively. For the sinusoidal rise, 
J = Jmsiti . (7rt/2r r ), while for the linear rise, J = Jiurt/Tr, where Jm is the maximal exchange couping during the 
pulse. 
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In Fig. |2| we plot the effective infidelity Aa z as a function of the pulse rise/fall time when the maximum exchange 
is set at 0.2meV. There are four sets (each with two pulse shapes) of data shown in the figure corresponding to four 
different values of Zeeman inhomogeneity, as indicated in the legend. At small pulse rise/fall time r, the infidelity 
Aa z grows approximately as a quadratic function of r (below ~ lOps) for both pulse shape. At longer t the growth 
saturates since the maximum value of Aa z is 2. In Fig. ^ we plot the effective infidelity Aa z as a function of the 
Zeeman inhomogeneity at a fixed exchange coupling J = 0.2meV and three different pulse rise times. Similar to Fig.|3 
the infidelity grows quadratically (for both pulse shapes) with inhomogeneity, then saturates when the inhomogeneity 
is in the same order of magnitude as the exchange coupling. The quadratic behavior in both Fig. [3] and Fig. 0] might 
be understood as the result of small parameter Taylor expansions. In Fig. [3 the infidelity Aa z is plotted as a function 
of the maximum exchange coupling J at a fixed inhomogeneity £7 = O.OlmeV and three different pulse rise times. 
Here for the two different pulse shapes the infidelity decreases with different power law as the exchange coupling J 
increases, essentially because the larger exchange coupling leads to an enhanced role played by the time-dependent 
function of the pulse shape. 

It is clear from Figs. to that imperfect exchange pulses (or non-square pulses) in the presence of Zeeman 
inhomogeneity reduce the effectiveness of the exchange interaction in performing swap operation. The key question 
now is whether one can still perform swap operations when square pulses are not available. The answer is affirmative. 
In analogy to the square pulse case shown in Fig. ^ one can turn on the most efficient exchange pulse as shown 
above, then let the system freely evolve (in the presence of Zeeman inhomogeneity) an optimal period of time, then 
turn on an exchange pulse again. Now each of the pulses is generally not an exact 7r pulse, but have to be calculated 
numerically to last for the most efficient duration. Similar to the square pulse case, such pulse sequence might have 
to be used more than once if the Zeeman inhomogeneity is too large while the exchange coupling strength is weak. 
However, in many realistic situations, some numerical corrections to the pulse sequence shown in Fig. ^ should be 
sufficient to produce a precise swap operation. In Tabled we summarize three different situations. Here we can see 



Pulse sequence 


1 


2 


3 


Zeeman inhomogeneity (meV) 


0.1 


0.1 


0.01 


pulse rise time (ps) 


5 


10 


5 


1st exchange coupling (meV) 


0.2 


0.2 


0.2 


2nd exchange coupling (meV) 


0.0965 


0.677 


0.000605 


corrections to the 2nd pulse (meV) 


+0.0465 


+0.672 


+0.000105 


fidelity 


0.999998 


0.999994 


0.9999999 



TABLE I: Numerically obtained pulse sequences for various exchange pulses to achieve a complete swap. Notice that in each 
of the cases the strength of the second exchange pulse increases, by 93%, 12.5 times, and 20%, respectively. 

that longer pulse rise time leads to significantly higher requirement for large exchange interaction. In the case of the 
first and second pulse sequences, if the pulses are square, the magnitude of the second exchange pulse needs to be 
S 2 / J — 0.05 meV. When the pulse rise time is 5 ps, this magnitude nearly doubles; while a 10 ps rise time leads to 
a large exchange coupling of 0.677 meV, more than ten times bigger than the square pulse case. If the maximum 
exchange is below 0.5 meV, a swap would require several pulses to perform (as in Eq. ©), so that the operation 
becomes less efficient. To further illustrate this point, we plot the strength of the second exchange pulse J' in the 
3-pulse swap sequence of Eq. (JSJ as a function of the pulse rise time in Fig. |SJ The unit for J' is its value in the 
square pulse case: 5 2 /Jq where S is the Zeeman inhomogeneity and Jo is the strength of the first exchange pulse 
(with the same pulse rise time for simplicity). It is quite clear from Fig. that the required strength for the second 
exchange pulse increases exponentially with the pulse rise time. For example, at 15 ps pulse rise time, the required 
J' is above 2.8 meV (as compared to 0.05 meV in the square pulse case), beyond the range available in the current 
state-of-the-art double dots. Multi-pulse sequences as in Eq. © would have to be invoked in such scenarios. 

In summary, we have assessed the effects of qubit inhomogeneity on the swap operations in a quantum dot quantum 
computer. We have shown that the exchange Hamiltonian becomes less efficient in performing swap when qubit 
inhomogeneity is present and the exchange is not turned on and off in the ideal square pulse shape. We have 
demonstrated ways to perform complete swap, at the expense of longer pulse sequences and numerically-searched 
pulse parameters |l8j . In exchange-based quantum computing schemes, entangling operations such as controlled- 
NOT are constructed from the square root of swap operation. With the extra complexity in the swap operation in 
the currently studied situation, it is quite natural to expect more complexity in controlled-NOT operations as well. 
Furthermore, since our calculation is performed in an effective two-level system, our results are applicable to single 
qubit operations, too. In essence, if |0) and |1) states are split energetically, attempts to perform direct rotations 
around axes other than z will be hindered if that interaction is turned on and off gradually instead of in a square 
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pulse profile. Indeed, the effects of such non-ideal pulses have been ex plor ed in solid-state systems like Cooper pair 
boxes 0, 0, |2(| and charge oscillations in double quantum dots pH \2j |23| , and will surely be encountered even 
more in the future experimental studies of various solid-state quantum computing architectures. 
This work is supported by ARDA and LPS. 
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FIG. 1: Square pulse sequences for a swap operation, (a) In the absence of Zeeman inhomogeneity, swap can be achieved by a 
7T pulse of the exchange interaction; (b) In the presence of Zeeman inhomogeneity 5, swap can be achieved by two n pulses of 
the exchange interaction with magnitudes J and J' = 8 2 /J, and a ir pulse for free evolution with S in between. 




FIG. 2: Sample trajectories of rotations on the Bloch sphere of the effective spin-1/2 system performed by single exchange 
pulses and projected onto the xy plane. 
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FIG. 3: Single pulse swap infidelity A<r 2 as a function of pulse rise time at a fixed exchange coupling J and four different Zeeman 
inhomogeneities. In this figure and the two following ones, the dashed lines without symbols represent data for trapezoidal 
pulse shape (linear rise/fall), while the solid lines with symbols are for a sinusoidal pulse rise/fall. 
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FIG. 4: Single pulse swap infidelity Aa z as a function of Zeeman inhomogeneity at a fixed exchange coupling J and three 
different pulse rise times. 



8 



10" 



< 



10" 



10" 



10" 



5y= 0.01 meV 



▼ ▼▼▼▼▼▼ 




time 10 ps 



rise time 10 ps 
30 ps 
50 ps 



10"' 
J (meV) 



10" 



FIG. 5: Single pulse swap infidelity Aa z as a function of the exchange coupling J at a fixed Zeeman inhomogeneity and three 
different pulse rise times. 




FIG. 6: Strength of the second exchange pulse in a complete swap pulse sequence as a function of the pulse rise time. The 
first exchange pulse strength is set at 0.2 meV while Zeeman inhomogeneity is 0.1 meV. 



